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Propagation of extragalactic ultra-high energy cosmic-ray nuclei : implications for 

the observed spectrum and composition 
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(— I ' In this paper, we study the propagation of cosmic-ray nuclei and protons. We emphasize the 

influence of the source composition on the expected spectrum and composition on earth as 
I ' well as on the phenomenology of the transition from Galactic to extragalactic cosmic-rays and 

O ' the interpretation of the ankle. We point out that the different source composition models 

cannot be distinguished on the sole basis of the cosmic-ray spectrum but that the the energy 
evolution of (Xmax) should remove this degeneracy. Finally, we compare the predictions for 
the different source composition models with the available data and discuss the implications 
of our results. 
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lO . 1 Introduction 

\ The cosmic-ray (CR) spectrum has been measured for half a century over 12 order of magnitude 

' in energy and 32 in Flux. This spectrum shows an extraordinary regularity as it can almost 

■ be fitted by a single power law between 10^^ and lO^'^ eV. One can however notice that this 

regularity is broken at least in two places : the so-called "knee" (a steepening around 3-5x10^^ 
and the "ankle" (a hardening around 3-5x10^^ eV). These features are of considerable 
^ ' interest for astrophysics. At the highest energies a suppression of the flux seems to be observed 

independently by two experiments (HiR es*^ an d Auger 131) above 4 — 5 x 10^^ eV probably 
corresponding to the expected GZK cut-offl23ISl_ However, at high energy, composition analyses 
are difficult to perform due to the large hadronic model dependance of air showers properties. 
As a result, a clear understanding of the highest energy cosmic-ray spectrum and its composition 
cannot yet be obtained. In particular, the question of the transition from Galactic (OCR) to 
extragalactic cosmic-rays (EGCR) is under intense debates. Moreover, a clear identification of 
the UHECR sources, crucial to constrain the acceleration mechanisms at play is not possible 
yet. In this context, propagation studies of extragalactic ultra-high energy cosmic-ray (UHECR) 
nuclei can be very useful to identify spectral and composition signatures of different source 
composition models that can be compared to the available data. 

In this proceeding, we study the propagation of cosmic-ray protons and nuclei. After present- 
ing the relevant photon backgrounds and interaction processes relevant at the highest energies, 
we calculate propagated spectra and the expected energy evolution of (Xmax) under different 
assumptions for the source composition and cosmological evolution of the source luminosities. 
We finally compare the model predictions with recent cosmic-ray data and discuss our results. 



2 Interactions of protons and nuclei with photon backgrounds 



Protons and atomic nuclei are usually considered as the most likely primary particles to explain 
the UHE cosmic-ray fluxes observed on Earth. Since the pioneer study of Puget, Stecker and Bre- 
dekamjP, many authors have considered the propagation of UHE nucl j^ ' ^™^ ' ^" ' ^^ ' ^^ ' ^'^ ' ^^ ' ^^ ' ^^ ' ^'^ ' ^^ ' ^^ ' ^" ' ^^ 

In the following sections we consider the interaction of protons and nuclei with the CMB and 
the infra-red, optical and ultraviolet backgrounds (hereafter we group these three backgrounds 
under IR/Opt/UV for short). To model the IR/Opt/UV backgrounds and their cosmological 
evolution, we use the latest estimate of^^H which is based on the earlier work of the authors 
updated with recent data on history of the star formation rate and the evolution of galaxy 
luminosity functions. We use IR/Opt/UV calculated at 26 different redshifts (Az = 0.2) between 
and 5. 

Protons and nuclei propagating in the extragalactic medium interact with CMB and IR/Opt/UV 
background photons. T hese in teractions produce features in the propagated UHECR spectrum 
such as the GZK cutoffESEl 

and their decay products generate the cosmogenic neutrino flux 
^23. In the case of protons the energy losses are dominated at low energy by adiabatic losses. 
Interactions with CMB photons become relevant at ~ 10^^ eV (at z=0) through the pair pro- 
duction process, dominant up to ~ 7 x 10^^ eV where the pion production, responsible for the 
GZK suppression, takes over. Interaction of protons with IR/Opt/UV photons are subdominant 
on the whole energy range. 

The interactions experienced by nuclei with photon backgrounds are different from the proton 
case. Pair production (for which we use the mass and charge scaling given in^ and adiabatic 
losses result in a decrease of the Lorentz factor of the UHE nucleus, whereas photodisintegration 
(also called photoerosion) processes lead to the ejection of one or several nucleons from the 
nucleus. Different photoerosion processes become dominant in the total interaction cross section 
at different energies^. The lowest energy disintegration process is the Giant Dipole Resonance 
(GDR) which results in the emission of one or two nucleons and a particles. The GDR process 
is the most relevant as it has the highest cross section with thresholds between 10 and 20 MeV 
for all nuclei. For nuclei with mass A > 9, we use the theoretically calculated GDR cross 
sections presented in^^S]^ which take into account all the individual reaction channels and are 
in better agreement with data than previous treatments. For nuclei with ^ < 9, we use the 
phenomenological fits to the data provided by^^l. Around 30 MeV in the nucleus rest frame and 
up to the photopion production threshold, the quasi-deuteron (QD) process becomes comparable 
to the GDR and dominates the total cross section at higher energies. The photopion production 
(or baryonic resonances (BR)) of nuclei becomes relevant above 150 MeV in the nuclei rest frame 
(e.g., ~ 5 X 10^^ eV in the lab frame for iron nuclei interacting with the CMB), and we use the 
parametrization given in^ where the cross section in this energy range is proportional to the 
mass of the nucleus (nuclear shadowing effects are expected to break this scaling above 1 GeV). 
The reference for this scaling is the deuteron photoabsorption cross section which is known in 
great detail. 

The contribution of the different photoerosion processes and the different backgrounds to 
the total mean free path for iron nuclei are displayed in Fig. [1^. The photoerosion is dominated 
by the GDR process through most of the Lorentz factor range. The baryonic resonances begin 
to dominate only above lO'^^'^ eV where the effect of the GDR starts to decrease. Fig. [T}) shows 
the contribution of pair production and photoerosion processes to the total attenuation length 
of iron nuclei. Photoerosion processes dominate through most of the energy range and the effect 
of pair production is small at low redshifts. Although the competition between pair production 
off the CMB and photoerosion processes with IR/Opt/UV photons depends on the redshift 
(e.g., at high redshifts pair production increases due to the stronger evolution of the CMB), the 
propagation of nuclei is mainly dominated by photoerosion processes. A comparison between 



the attenuation lengths of different species is displayed in Fig. [T];. The figure shows what is 
known since^^, that the attenuation length of low mass nuclei are smaller than that of protons 
and heavy nuclei above 10^^ eV and, as a consequence, light nuclei should not contribute as 
significantly at the high energy end of the spectrum. Furthermore, iron nuclei have larger or 
similar attenuation lengths to protons up to ~ 3 x 10^'' eV. However, the energy loss processes are 
different for protons and nuclei and the sole comparison of attenuation lengths can be misleading 
and is not straightforward to interpret : most of the energy losses of nuclei result in nucleon 
ejection, thus, unlike protons, a given nucleus does not remain on "the same attenuation length 
curve" during its propagation. 
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Figure 1: Left: Evolution of the iron nucleus mean free path as a function of the Lorentz factor for the different 
photoerosion processes and interactions with the CMB and IR/Opt/UV photons at z = 0. Center: Evolution of 
the attenuation length at z = 0. The contribution of pair production and photoerosion processes off the CMB 
and IR/Opt/UV photons are separated. Right: Comparison of the energy evolution of the attenuation length of 

different nuclei at z = 0. 



3 Propagated spectra 

For the following calculation we use a M onte- Carlo code to propagate nuclei from the source 
to Earth as described in detail in^^ZISS!^ "\Yg consider the classical pure proton scenarios and 
the extragalactic mixed composition models. For our generic mixed nuclei case (see^^^ for more 
details), we assume that the EGCR source composition matches that of the GCRs observed at 
lower energies, and that the maximum energy achieved by nuclei of species i in EGCR sources 
scales with their charge Zj, i.e. E^^^^^i = ZiE^^^i^H), as expected if the acceleration mechanism 
is controlled by magnetic confinement and limited by particle escape. In the following, we 
assume a power-law source spectrum with spectral index (3 and set the maximum proton energy 
to lO^"'^ eV, unless otherwise specified, referring to"^^ for a discussion on the influence of Emax- 
We find that the observed UHECR spectra are best fitted with spectral indices between 2.1 
and 2.3, which corresponds to a proton dominated composition with significant fractions of He 
and CNO, and a lower fraction of heavier nuclei^^Sl. However, another important ingredient of 
EGCR models is the time evolution of the power and/or number density of sources. Indeed, 
the link between the spectrum of the sources and the observed one (and thus the determination 
of the "best fit spectral indices") depends strongly on the assumed redshift evolution of the 
sources. Here, we consider three different source evolution models. The first one corresponds 
to no evolution at all - hereafter referred to as the uniform source distribution model. In the 
so-called "SFR model", we assume that the EGCR injection power is proportional to the star 
formation rate, which correspond a to redshift evolution in (1 -|- z)^ for z < 1.3 and a constant 
injection rate for 1.3 < z < 6 (with a sharp cutoff at z = 6). Finally, we consider a stronger 
source evolution favoured by the recent infra-red survey of the Spitzer telescope. In this so-called 
"strong evolution model", we assume a injection rate proportional to (1 -|- z)'* for z < 1 and a 
constant rate for 1 < z < 6, followed by a sharp cut-off (see^^El for more details and references 




Figure 2: Propagated spectra, E^^{E), for pure proton models (left) and mixed composition models (right) 
compared with HiRes monocular data. Different source evolution models are indicated by the labels. The 
corresponding galactic components are inferred from the overall spectrum by subtracting the EGCR component, 
in the case of the uniform and SFR source evolution models. 



on the sources evolution models). 

In the case of pure proton EGCR sources, the best fit /3 = 2.6 if one assumes a uniform 
distribution of sources (no evolution) , while it goes down to 2.5 in the case of an SFR-like evolu- 
tion, and 2.4 in the strong evolution case (see Fig. Hi and^^^^. As shown in previous works 
l ^^ l "^^ !, the propagated proton spectrum and the concave shape known as the "pair production 
dip" (with a minimum around 10^^ '' eV on Fig. [2^, for the uniform source model) are only 
mildly dependent on the source evolution hypothesis. However, the energy where this e''~-e~ dip 
begins depends on the relative weight of energy losses related to pair production, which domi- 
nate at high ener gv, an d energy losses associated with the universal expansion, which dominate 
at low energy more details). Therefore, the beginning of the dip depends on 

the redshift evolution of the source density (or power): the transition between the two energy 
loss processes occurs at a lower energy in the SFR and strong evolution cases (see Fig. [2^). 
In these cases the extragalactic component can account for the whole CR flux down to much 
lower energies (~ 410^^ eV), which correlatively allows/requires the GCR component to cut 
at relatively low energies, notably lower than the confinement limit of charged nuclei in the 
Galaxy. Thus, in the pure proton case, the energy Eend at which the GCR/EGCR transition 
ends (i.e, above which cosmic-rays are purely extragalactic) depends on the source evolution 
scenario. The highest value of Send is obtained in the case of a uniform source distribution, 
around 1-1.5 10^^ eV. This energy range is significantly lower than in the case of mixed compo- 
sition scenarios (see below) - a distinctive feature that can be used to discriminate between the 
models, using composition analyses. 

The propagated EGCR spectra obtained with a mixed source composition are shown in Fig. 
[2}3. The best fit of the high-energy data is obtained in these cases for significantly smaller spectral 
indices, i.e., harder source spectra: P ~ 2.3 in the uniform case, going down to 2.2 for SFR- 
like source evolution and 2.1 for the strong evolution model. In all these mixe d corn position 
cases, the end of the GCR/EGCR transition roughly coincides with the ankle ^2Zl3Qj Above 
j^gis.s i-j-^g predicted spectrum is quite insensitive to the source distribution, as can be seen 
in Fig. [2)3. It is also important to note that the mixed composition models do not imply/require 
any definite value of the highest energy of cosmic rays in the Galactic component, as long as 
GCRs represent a sufficiently small fraction of the total spectrum around Eankie not to influence 
the overall spectrum and composition. Therefore, the Galactic component does not necessarily 
vanish above -Bankie) nor is it required that cosmic rays be accelerated above -Eankie at all. At 
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Figure 3: Left: (Xmax) evolution for an extragalactic mixed composition and uniform source evolution model for 
three different hadronic models (see labels). Right: Predicted spectrum for a mixed extragalactic composition 

above 10^^ eV decomposed in its elemental components. 



energies below the ankle, the inferred fraction of GCRs depends on the source evolution model, 
just as in the pure proton case (see'^for more details). We stress that the energy at which the 
Galactic and extragalactic components have an equal contribution to the CR flux lies between 
~ 5 10^"^ eV and ~ 10^^ eV . Note also that in our mixed composition models the cosmic rays 
can be dominated by light nuclei at energies below 10^® eV, which is a major difference with the 
GCR/EGCR transition scenario studied in^^. 



4 The shape of {X^,^,)(E) 

Having identified the value of (3 that provides the best fit of the data in each scenario, and 
obtained the corresponding fractions of GCRs and EGCRs at all energies, we can now deduce 
the evolution of the cosmic-ray composition as a function of energy and predict the values of the 
associated observables. The propagated EGCR composition in each case is a direct output of 
our computations, and we assume that the Galactic component is essentially made of Fe nuclei 
above 10^^'^ eV (relaxing this assumption would slightly fiatten the (Xmax) evolution in the 
transition region). From the relative abundance of all elements at a given energy, we derive the 
average value of the atmospheric depth (in g/cm^) at which the maximum shower development 
is reached, (Xmax)) using Monte-Carlo shower development simulations'^^ 

In the pure proton case, t he inte rpretation of the evolution of (Xmax) with energy is straight- 
forward (figures are shown inEGEHl^^ ^he transition from Galactic iron nuclei to extragalactic 
protons being quite narrow (i.e., it occurs over a small energy range, in a decade or even half a 
decade), the evolution of (Xmax) with energy is very steep and then gets flatter when the tran- 
sition is over and the composition does not change anymore, all EGCRs being merely protons. 
The point where the (^max) evolution can be observed to break simply indicates the energy £^end) 
corresponding to the end of the transition. Characteristically, an early break in the elongation 
rate at ~ 4 10^^ eV is expected in the strong and SFR source evolution models, whereas the 
break is found around 1-1.5 10^^ eV for a uniform source distribution. No break is expected at 
the ankle. Indeed, the ankle is consistently interpreted in pure proton models as the signature 
of the interactions between EGCR protons and CMB photons producing e"'"e~ pairs. Obviously, 
the resulting "pair production dip" would not be visible if the EGCR component did not consist 
almost exclusively of protons. Quantitatively, nuclei he avier than H cannot contaminate the 
EGCR component at a higher level than ~ 15 opZESET] 




Figure 4: Left: (Xmax) evolution for an extragalactic mixed composition above 10^^ eV, for the SFR and uniform 
source evolution compared to cosmic-ray data (see legend). Right: Predicted spectrum for a mixed extragalactic 
composition {Emax = Z x 10^^ eV) decomposed in its elemental components and compared to Auger data. 



The case of mixed composition models is illustrated in Fig. [3^. The evolution of (Xmax) is 
relatively steep in the transition region, below Eankie, because the composition evolves rapidly 
from the dominantly heavy Galactic component to the light extragalactic mixed composition. 
However, the evolution is significantly slower than in the case of pure proton models, because 
the transition is wider and the cosmic-ray composition does not turn directly into protons only. 
As can be seen on Fig. [3^, an intermediate stage appears, which may be called the mixed- 
composition regime, where a break in the evolution of (Xmax) around Eankie is followed by a 
flattening up to ~ 10^^ eV, reflecting the fact that the (propagated) EGCR composition does 
not change much in this energy range. This is because among the different EGCR nuclei, 
only He nuclei interact strongly with infrared photons at these energies. Between Eankie and 
~ 10^^ eV, the evolution of (Xmax) is actually compatible with what is expected from a constant 
composition. Then around 10^^ eV, the relative abundance of nuclei heavier than protons starts 
to decrease significantly as a result of photo-disintegration processes: the CNO component 
starts interacting with the infrared background and the CMB photons eventually cause the He 
component to drop off completely (see Fig. [Sb)- The evolution of (Xmax) therefore steepens 
again, accompanying the progressive evolution towards an almost pure proton composition as 
each type of nuclei reaches its effective (mass dependent) photo-disintegration threshold. Even 
though slight differences may be expected from one model to the other, the above evolution of 
(-'^max )(-£') in a three steps process is a characteristic prediction of mixed-composition models, 
or generically of any type of EGCR sources allowing for the acceleration of a significant fraction 
of nuclei heavier than He. 



5 Comparison w^ith data and discussion 

The high-energy cosmic ray spectrum can be satisfactorily accounted for within either the pure 
proton or the mixed composition models (and many other source composition models, see'^^^. 
However, we have shown that the corresponding phenomenology of the GCR/EGCR transition is 
very different in each case, which results in distinct shape of (Xmax) as a function of energy. The 
currently available data do not allow one to draw definitive conclusions yet. However, we argued 
in^S^that the predictions of the mixed-composition models appear to be in better agreement with 
the current data from fluorescence detectors. In particular, a good agreement is found with Fly's 
Eye results above 10^'^'^ eV^''^-. Concerning the slope of the (Xmax) evolution in the transition 
region (i.e., below the ankle), mixed-composition models typically predict values for the slope of 



the (Xmax) evolution compatible with what is observed. Furthermore, both the predicted break 
at the ankle and the steepening above 10^^ eV are compatible with the HiRes Stereo data^^. 
Fig. shows the comparison between the (Xmax) evolution for a mixed composition in the SFR 
evolution case and the data of HiRes Stereo, HiRes- Mia, Fly's Eye (rescaled by 13 gcm~'^, as 
suggested inS31) Qjy(\ AugerHS. As can be seen. Fly's Eye and Stereo HiRes data are consistent 
with the predicted break in the (Xmax) evolution between 3 and 4 EeV, which also corresponds 
to the energy of the ankle reported by both experiments. Note that HiRes-Mia results at lower 
energy are also compatible with pure proton models, as well as with mixed-composition models 
except at one point around 5 10^'^ eV. The data best agree with the absolute scale of (^max) 
computed with the QGSJet-H model. However, it is important to note that the results obtained 
with QGSJetOl show exactly the same features shifted downwards by ~20 gcm~'^, and are still 
well within the systematic uncertainties of the different experiments. This illustrates once more 
that the choice of the hadronic model is not critical in the present discussion's^. 

When comparing the mixed composition model predictions with the recent data of the Pierre 
Auger observatory, the agreement appears less good. Auger data are compatible with a break 
in the (Xmax) evolution in an energy range close to the ankle (which is, as mentioned before, 
difficult to handle for a pure proton model) and the overall shape of the evolution below 10^^ eV, 
although flatter, is compatible with the mixed model predictions. Above ~ 10^^ eV, however, 
the composition seems to be getting heavier where the mixed model would predict a lightening. 
This trend, if confirmed by a larger statistics, would then represent a major incompatibility 
with our "generic" mixed composition model. Though the reality of this trend is still quite 
speculative due to the low statistics accumulated at the highest energies, it is interesting to 
investigate what could cause the composition to get heavier at the highest energies. 

In the case of our generic mixed composition, one actually expects that the composition 
could get heavier above 5 x 10^^ eV if heavy nuclei are accelerated at the highest energies. 
Indeed, above energy, the proton flux decreases sharply due to photopion interactions with 
CMB photons. Between 5 x 10^^ eV and ~ 2.5 x 10^° eV, heavy nuclei (Fe) only interact 
with the less dense far-IR photons and the decrease of their flux is then slower (see Fig. [3b) 
than for protons, implying an increase of their relative abundance. This trend ends when the 
heavy component disappears due to interaction with CMB photons (see^^Sl for more details). 
However, the energy range where this feature is expected does not seem to be compatible with 
Auger data. Invoking a heavy dominated or possibly pure iron composition at the sources would 
presumably not solve the problem either. The composition analyses at lower energy seem to 
favor a transition from heavy galactic to light extragalactic cosmic-rays which does not look 
compatible with the expectations of an extragalactic heavy source composition (see discussion 

on 

of the pure iron sources in"^^). Let us note that the latter point is somehow alleviated but 
remains true if one assumes a strong evolution of source luminosity with redshift for which 
spectral indices of 2.0-2.1 would be required to fit the observed spectra. 

In this context, a possible explanation for this trend would be to infer that most of the 
sources are not able to accelerate protons up to the highest energies and that only heavy nuclei 
are accelerated above 10^^ eV (in the case of a maximum energy scaling with the charge of the 
nucleus, this kind of scenario is proposed for instance in^Sl). As we discussed in^SQl^ fo^ E^ax 
proton solutions do not work very well with our usual mixed composition hypothesis. Some 
tuning of the composition is then necessary for this type of scenario to be compatible with the 
data. However, good fits of the data can be obtained with a moderate increase of the overall 
abundance of heavy elements at the source, by a factor of three or so with respect to low-energy 
Galactic cosmic rays. This is illustrated in Fig. db where expected spectra are displayed and 
compared with data assuming a mixed composition, /3 = 2.0 , Emax = Z x 10^^ eV and ~ 30% 
of Fe nuclei at the sources. One can see that the agreement with data is reasonable and that the 
composition, proton dominated at low energy, becomes gradually heavier and very dominated 



by iron above 5 10 eV. Such a scenario would have imphcations on the cosmogenic neutrino 
flux that should be extremely low above 10^'^ eV. Neutrinos and photons fluxes from interactions 
during the acceleration at the sources would as well be presumably very low (see discussion in 
To conclude, let us note that, at the current level of statistics, it seems difficult to argue that 
the recent anisotropy claim by the Pierre Auger collaborationSSl disfavors heavy compositions at 
the highest energie (see for instance^^^. However, a significant small scale clustering, if observed 
in the future, would challenge this type of scenario. 

Although the available cosmic-ray data already allow to put encouraging constrains on the 
different models we presented, a clear picture of the GCR to EGCR transition, the composi- 
tion and the origin of the UHECR is still difficult to draw. Future data of the Pierre Auger 
Observatory, its low energy extension as well as the future very large aperture projects such as 
Auger 

NorthlSni and JEM-EUSoEI should allow a clearer view on these questions in the next 
few years. 
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